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Abstract 

We study CP asymmetries in two-body decays of bottom squarks into charginos 
and tops. These asymmetries probe the SUSY CP phases of the sbottom and 
the chargino sector in the Minimal Supersymmetric Standard Model. We iden- 
tify the MSSM parameter space where the CP asymmetries are sizeable, and 
analyze the feasibility of their observation at the LHC. As a result, potentially 
detectable CP asymmetries in sbottom decays are found, which motivates fur- 
ther detailed experimental studies for probing the SUSY CP phases. 
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1 Introduction 



Supersymmetry (SUSY) P is a well motivated theory to extend the Standard 
Model (SM) of particle physics. SUSY models are not only favored by gauge cou- 
pling unification and naturalness considerations, but are also attractive from the 
cosmological point of view. For instance the lightest SUSY particle (LSP) is a good 
dark matter candidate if it is stable, massive and weakly interacting [2l|3]. SUSY 
models can also provide new sources of CP violation In the Minimal Super- 
symmetric Standard Model (MSSM) [1], the complex parameters are conventionally 
chosen to be the Higgsino mass parameter /i, the U(l) and SU(3) gaugino mass 
parameters Mi and M3, respectively, and the trilinear scalar coupling parameters 
Af oi the third generation sfermions (/ = 6, t, r), 

/i = |/^|e^'^^ Ml = |Mi|e^'^S M3 = {M^le''^', Af = \Af\e'^^f . (1) 

These phases contribute to the electric dipole moments (EDMs), in particular to 
those of Thallium [5], Mercury |6], the neutron [^, and the deuteron |8], which can 
be beyond their current experimental upper bounds [SHS]. The experimental limits 
generally restrict the CP phases to be smaller than vr/lO, in particular the phase 
[TOtfTT] - However, the extent to which the EDMs can constrain the SUSY phases 
strongly depends on the considered model and its parameters [TOl - fTG] . 

As shown for example in Ref [16], due to cancellations among different con- 
tributions to the EDMs, only isolated points in the CP phase space can give large 
CP-violating signals at the LHC. It is important to search for these signals, since the 
cancellations could be a consequence of a deeper model that correlates the phases. 
In addition, the existing EDM bounds can also be fulfilled by including lepton flavor 
violating couplings in the slepton sector [T3j. This is important when considering for 
example SUSY Seesaw models, where CP violation in the slepton sector is connected 
to the neutrino sector and Leptogenesis pTlfTB] . 

Thus measurements of SUSY CP observables outside the low energy EDM sector 
are necessary to independently determine or constrain the phases. In particular, the 
phases of the trilinear scalar coupling parameters Af have a significant impact on 
the MSSM Higgs sector [19]. Loop effects, dominantly mediated by third generation 
squarks, can generate large CP-violating scalar-pseudoscalar transitions among the 
neutral Higgs bosons [SUIET]. As a result, the lightest Higgs boson with a mass of 
order 10 GeV or 45 GeV [22] cannot be excluded by measurements at LEP [23]. 
The fundamental properties and the phenomenology of CP-violating neutral Higgs 
boson mixings have been investigated in detail in the literature [2^125]. 

The phases can also drastically change other SUSY particle masses, their cross 
sections, branching ratios [261429] . and longitudinal polarizations of final fermions [30] . 
Although such CP-even observables can be very sensitive to the CP phases (the ob- 
servables can change by an order of magnitude and more), CP-odd (T-odd) observ- 
ables have to be measured for a direct evidence of CP violation. CP-odd observables 
are, for example, rate asymmetries of cross sections, distributions, and partial decay 
widths [31]. However, these observables require the presence of absorptive phases. 
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e.g. from loops. Thus they usually do not exceed the size of 10%, unless they are 
resonantly enhanced [2I1I2S1IS21E3]- 

Larger CP asymmetries in particle decay chains can be obtained with triple 
products of final particle momenta [M]. They already appear at tree level due to 
spin correlations. Triple product asymmetries have been intensively studied in the 
production and decay of neutralinos [351439] and charginos [391 - H2] at the ILC [13], 
also using transversely polarized beams [H]. At the LHC [15], triple product asym- 
metries have been studied for the decays of neutralinos [STlllSlllZ] , stops [THlllHlllH] , 
and sbottoms [50]. For recent reviews, see Ref. \5T\ . 

In this paper, we thus study CP asymmetries in two-body decays of a sbottom, 
b^-^t + x~; m = l,2; z = l,2; (2) 
followed by the subsequent two-body decay of the chargino, 

X7^^i+^h ^ = e,/i, (3) 
with the invisible sneutrino decay i/| — Xi ^i, and the subsequent top decay 

t^b + W-^; W+^ue + it; i = e,i^, (4) 

see Figure [T] for a schematic picture of the entire sbottom decay chain. The CP- 
sensitive spin-spin correlations of the sbottom decay allow us to probe the phase of 
the coupling parameter Af,, and the phase of the higgsiono mass parameter fi. The 
phases of the trilinear scalar coupling parameters of the third generation sfermions 
are rather unconstrained by the EDMs [TOllTS] . Therefore it is appealing to study CP 
asymmetries in squark decays at high energy colliders like the LHC [15] or ILC [13] . 
The third generation sfermions also have a rich phenomenology due to a sizeable 
mixing of left and right states. 

CP asymmetries based on triple products have already been studied in two-body 
decays of sbottoms [50], however in their rest frame only. At colliders like the LHC, 
the particles are highly boosted, which will generally reduce the triple product asym- 
metries [T6ll46p47[l^ . We thus will also include the sbottom production at the LHC. 
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In addition, we will also study asymmetries which base on epsilon products. Those 
asymmetries are boost invariant, and thus provide the largest possible asymmetries. 
We calculate the amplitude squared for the entire sbottom decay in the spin-density 
matrix formalism [52]. The compact form of the amplitude squared allows us to 
identify the optimal CP observables for sbottom decays, an asymmetry, that has 
not been studied in Ref. [50] . 

In Section [SJ we identify the CP-sensitive parts in the amplitude squared, de- 
fine the CP asymmetries in bottom squark decays, and discuss their dependence 
on the complex sbottom-top-chargino couplings. We analyze their MSSM parame- 
ter dependence in a SUSY benchmark scenario. In Section [3l we discuss sbottom 
production and their boost distribution at the LHC. We give lower bounds on the 
required LHC luminosities to observe the asymmetries over their statistical fluctu- 
ations. We summarize and give our conclusions in Section |H In the Appendix, we 
review sbottom mixing with complex parameters, give the phase space, and calculate 
the sbottom decay amplitudes in the spin-density matrix formalism. 



2 CP asymmetries in bottom squark decays 

In this Section, we identify the CP-sensitive parts in the amplitude squared of the 
entire two-body decay chain of the bottom squark, see Eqs. ([2])-(lll) and Figure [TJ 
In order to probe these parts, we define CP asymmetries of epsilon and triple prod- 
ucts of the particle momenta. Explicit expressions for the squared amplitude, La- 
grangians, couplings, and phase-space elements are summarized in the Appendix. 



2.1 T-odd products 

The amplitude squared |Tp for the sbottom decay chain, see Figure [H can be 
decomposed into contributions from the top spin correlations, the chargino spin 
correlations, the top-chargino spin-spin correlations, and an unpolarized part, see 
Eq. flE.SOp . Since we have a two-body decay of a scalar particle, a CP-sensitive part 
can only originate from the spin-spin correlations. It is given by the last summand 
in Eq. flE.55j) . which is proportional to 



|T|2 D Im{4i(fc^)*}[p^,pt,p^,,p^2], m,i = 1,2. (5) 

The left and right couplings and respectively, are defined through the bm- 
t-xf Lagrangian |18], 



^tl^+=9t{lLPR + kl,PL)xtK + ^-C: (6) 

see Appendix O These couplings depend on the mixing in the sbottom and chargino 
sector, and thus on the CP phases (pA^ and 0^. The imaginary part of the coupling 
product, Im{/^j(/i;^j)*}, in Eq. (|5]) is multiplied by a T-odd epsilon product for 
which we use the short hand notation 

£ = [Ph Pt^ Ph ' Pi2] = ^t^ua^ Pi Pt Ph pi ' C^) 
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with the convention £0123 = 1- Since each of the spatial components of the four- 
momenta changes sign under a naive time transformation, t — ?■ —t, this product is 
T-odd. Due to CPT invariance, T-odd products are related to CP-odd observables. 

2.2 T-odd asymmetries 

The task is to define an observable, that projects out the CP-sensitive part of the 
spin-spin correlation term from the amplitude squared. This can be achieved by 
defining for the T-odd product S, Eq. ([7]), the T-odd asymmetry of the partial 
sbottom decay width F [TB] , 



with the amplitude squared |Tp, and the Lorentz invariant phase-space element 
(iLips, such that J |Tp(iLips/(2m^) = F. The T-odd asymmetry is also CP-odd, if 
absorptive phases (from higher order final-state interactions or finite-width effects) 
can be neglected [M]. 

In general, largest asymmetries are obtained by using the epsilon product S, see 
Eqs. dS]) and ([7]), that matches the kinematic dependence of the CP-sensitive terms in 
the amplitude squared. In the literature, this technique is sometimes referred to opti- 
mal observables [53]. Other possible combination^ of momenta, £ = [Pi,,Pb,Pii,Pi2]^ 
or £ = [pi,Pt,Pb,Pii], lead to smaller asymmetries, see Ref. [50], and Fig. |6l 

Triple products of three spatial momenta can also be used to define asymme- 
tries [3ll|50]. In the sbottom rest frame, = (m^, 0), the epsilon product is 



That triple product will give the largest asymmetries in the sbottom rest frame, and 
the other combinations of momenta for T lead to smaller asymmetries. 

Note that the asymmetries of an epsilon product £ are by construction Lorentz in- 
variant whereas those constructed with a triple product T are not [TBll^ . The triple 
product asymmetries will therefore depend on the sbottom boost, f3i = \pi\/Ei, and 
are generally reduced if not evaluated in the sbottom rest frame. We will discuss 
the impact of the sbottom boost on the asymmetries at the LHC in Section |3l 



Note that momenta from both the decay products of the chargino and the top have to be 
inchided to obtain non- vanishing asymmetries [SD] . Otherwise the CP-sensitive top-chargino spin- 
spin correlations are lost. If for example momenta from the top decay are not taken into account, 
only CP asymmetries from the chargino decay can be obtained. In addition, a three-body decay 
is required, (or a two-body decay via an on-shell W boson [3D]) to probe the phases 0^, 0i of the 
chargino/neutralino system. The asymmetries are then of the order of 10%, which is typical for 
chargino three-body decays [32], and also neutralino three-body decays [5S1I551HBII471H5] . 



A 



r{£ > 0) - r{£ < 0) _ /Sign[g]|Tpc/Lips 
T{£ > 0) + F(£' < 0) ~ / |T|2rfLips 




(9) 
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Table 1: MSSM scenario. The mass dimension parameters are given in GeV. 
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2.3 Parameter and phase dependence 

In order to analyze the phase and MSSM parameter dependence of the asymmetry 
A, Eq. dH]), we insert the explicit form of the amplitude squared in the spin-density 
formalism. As shown in Appendix |Fl we obtain 

/ Sign(g)(pb ■ p^J [Pi,Pt,Ph,P£2] dUps 
{Pxt -Ph) I (Pt ■ Pi2 ) (Pb ■ Put ) c^Lips 

with ip ± ■ ptj = (m^± — m? )/2, and the coupling function 



Im{C(/cL)*} 
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1(1/^ |2-U|P |2) ■RpI/'' (k'' Y\ 



The asymmetry can thus be separated into a kinematical part and the effective 
coupling factor 77. That factor is approximately independent of the particle masses, 
but governs the main dependence on the CP phases, and on the parameters of the 
sbottom-top-chargino couplings /^j, fc^j. To qualitatively understand this depen- 
dence, we expand 



Im{/L(A;L)*} = c„rjm{f/* V;;} -\d^YtYbM2e-,)\m{U:,V:^e-''t'-^}, (12) 



with the Yukawa couplings Yt and Y^, see Eq. f lC.14p . the short hand notation 



Ci = cos^6'^, C2 = sin^6'^, di = 1, ^2 = —1, the sbottom mixing angle Eq. (lA.Sp . 
and the CP phase 0^ = arg[Ab — /i*tan/3], Eq. (IA.4p . of the sbottom system. The 
imaginary part of the products of the sbottom-top-chargino couplings Im{/J^j(fc^j)*} 
is sensitive to the phases 0^ and (pA^, and can be large due to the Yukawa couplings 
Ft, Ift. For 0^ = or tt, the chargino diagonalization matrices t/, V are real, and [50] 



oc sin(2%)sin(0^), (13) 

In particular, for 0^ = or vr, Im{/^-(A;|^J*} shows a sin(0yij behaviour, and 
is maximal at sin(0^j^) ^ 7r/2,37r/2. However, the maxima of the asymmetries 
correspond to the maxima of 77, which are shiftet away from (pA^ = vr/2, 37r/2. This 
is due to the influence of the denominator of rj, see Eq. flF.65p . which has a cosine-like 



dependence on (pA^y and 77 will be approximately maximal for f» l/c^^. 
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Table 2: SUSY particle masses for the benchmark scenario of Table [H = e, /i). 



= 400 GeV 


m^^ = 424 GeV 


m + = 158 GeV 


m + = 301 GeV 

A2 


m^o = 108 GeV 


m^o = 172 GeV 


= 206 GeV 


m^o = 302 GeV 


rrii}^ = 137 GeV 


TUf^ = 152 GeV 
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m^^ = 157 GeV 



2.4 Numerical analysis 

To quantitatively study the asymmetry A, Eq. ([H]), we analyze the behaviour of 
the coupling factor t] in more detail, and define a benchmark scenario in Table [TJ 
We give the relevant resulting SUSY masses in Table [2j We fix the soft-breaking 
parameters in the slepton sector by Af| = M| = 150 GeV for i = e, fi, r, to enable 

the subsequent chargino decay xf — )■ £f z>^*\ We take stau mixing into account, and 
fix the trilinear scalar coupling parameter Aj. = 250 GeV. Since its phase does not 
contribute to the CP asymmetry, we set it to (pA-r = 0, as well as (pi = 0, which is 
the phase of the U(l) gaugino mass parameter Mi. In order to reduce the number 
of parameters, we further use the GUT inspire 4l| relation | Mi | = 5/3 Ms tan^ 6'^. 

We will study the decay of the lighter sbottom into the lightest chargino, bi — 
txi, which gives the dominant contribution to the asymmetries. For our reference 
scenario as defined in Table [H the corresponding branching ratios are BR(&i — 
txf) = 19% and BR(x^ — e^Ue) = 33%. Other decay channels yield much 
smaller asymmetries for the scenario as defined in Table [T] Further, we will fo- 
cus on the largest (Lorentz invariant) asymmetry for the optimal T-odd product 
£ = [pi,Pt,Pei,Pe2]^ see Eq. ([7]). Note that in the sbottom rest frame, this asym- 
metry is equivalent to the triple product asymmetry with T = pt ■ {pe^ x pe^), see 
Eq. (|9]). We use the short hand notation A{tiii2) for the asymmetry, to indicate 
the momenta used for the triple product. 

In Figure El we show the dependence of the coupling factor t], Eq. f lTTj) . and its 
corresponding asymmetry A{tiii2), Eq. f lTU]) . on the two CP phases (pA^ and 0^, for 
the SUSY scenario of Tabled! We clearly see that the maxima of A{tiii2) ~ ±40% 
are not necessarily obtained for maximal CP phases, 0^^,^ = 7r/2, 37r/2. The reason 
is that the phase dependence of A is almost governed by the coupling factor r], shown 
in the left panel of Figure |2j Thus the asymmetry can be sizeable even for small 
values of the phases, favored by EDM constraints, which in particular constrain 
(p^. For example, the asymmetry has a maximum of A{tiii2) ~ ±40% at small 
(p^ ~ ±0.27r, for (pA^ = 0. The positions of the maxima will remain when including 
the effects of the sbottom boost, as we will discuss in Section l3m Figure [2] motivates 
the choice (pAi, = 0.27r, 0^ = 0, in our benchmark scenario. Table [T] 



"'Note that this choice also can significantly constrain the neutralino sector [57) . 
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Figure 2: Phase dependence of the CP-odd couphng factor rj, Eq. flTT]) . for sbottom 
decay bi — )■ txi (left), and the corresponding CP asymmetry A{tiii2), Eq. ([HD, 
in percent (right), for the subsequent two-body decay chain Xi ^ ~^ 

hvii2i see Figure [H in the hi rest frame. The SUSY parameters are given in Table [H 
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Figure 3: Contour lines in the yU-M2 plane of the CP-odd coupling factor r], Eq. (fTTj) . 
for sbottom decay bi txi (left), and the corresponding CP asymmetry A{tiii2), 
Eq. ([8]), in percent (right), for the subsequent two-body decay chain Xi ~^ 
and t biyii2, see Figure [H in the sbottom rest frame. The SUSY parameters 
are given in Table [TJ The area above the contour lines is kinematically forbidden 
by m^^ < m^± + nit, the area below the contour lines of the asymmetry (right) is 
forbidden by m ± < rriy^. The area above the dashed line is excluded by mf^ < m^o. 
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Mo [GeV] tan/3 

Figure 4: Contour lines of the CP asymmetry A{tlil2)i Eq- dH]), in percent for 
sbottom decay hi — )■ txi , followed by the subsequent two-body decay chain Xi ~^ 
z>|, and t — > 6 z/^ in the plane of the soft breaking parameters Mg, (left), 
and in the tan/3-|y4fe| plane (right). The SUSY parameters are given in Table [TJ 

In Figure ini we present the dependence of the asymmetry (left) and the coupling 
factor rj (right) on the chargino mixing, which is mainly determined by the higgsino 
and gaugino parameters /i and M2, respectively. We find large values of the coupling 
factor rj and the asymmetry for mixed and higgsino-like charginos. The kinematical 
boundaries for the sbottom decay bi tx^ , and the chargino decay x^ — )■ z>| are 
indicated in Figure [31 

The size of the asymmetry A, Eq. fITU]) . strongly depends on the sbottom mixing, 
see the mixing matrix Eq. (lA.ip . in Appendix |Al The mixing is determined by the 
soft-breaking parameters Mq, Mfj, on the diagonal, and mh{Ah — /i*tan/3) on the 
off-diagonal entries of the sbottom mixing matrix. In Figure H] (left), we see that 
the asymmetry is maximal in the region of the level crossing Mq ^ M^, where 

also lm{l\^{kl^)*} is maximal, see Eq. f[T^ . However, although the imaginary part 
quickly drops for increasing Mq > Mfj, the asymmetry is still sizeable in that region, 

since the couplings still fulfill {l^il ~ \k\i\. 

In Figure m (right), we show in the tan/3-|y4fc| plane contour lines of the asym- 
metry, which peaks at \Af,\ ^ tan/3, where ~ l^iil- Whereas the imaginary 
part of the couplings, lm{l\^{k\i)*} , steadily grows with increasing \Ai,\ > |/i|tan/3, 
see Eq. (1121) . the asymmetry is reduced in that region. This is since the sbottom 
width T(bi txi) increases, which enters in the denominator of the asymmetry, 
see Eq. ([8]). Thus from Figure HI we can observe that the contribution from the 
CP-even denominator of 77, Eq. (fTT]l . to the asymmetry has an important impact on 
its parameter dependence. 
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Figure 5: Total sbottom pair production cross section a(pp — )■ bb*) at the LHC 
as a function of the sbottom mass (left). Normalized sbottom pair production 
distribution with respect to the sbottom boost factor (right). The leading 

order cross sections have been calculated at y/s^ = 14 TeV using MadGraph [55] . 



3 CP asymmetries at the LHC 

The production of sbottom pairs at the LHC [H2] 

p + p-^K, + b*^, m = l,2, (14) 

dominantly proceeds via gluon fusion. As a result, the leading order cross section 
a{pp — )■ &m&m) is independent of any other SUSY model parameters than the sbot- 
tom mass. The strong dependence can be seen in Figure [5] (left), where the cross 
section drops by six orders of magnitude with an increase of the squark mass from 
0.2 TeV to 2 TeV. The produced sbottoms have a distinct distribution in their boost 

Pi = (15) 

along the direction of their momenta. In Figure [5] (right), we show the normalized 
boost distribution for several values for the sbottom mass. Typically, light sbottoms 
are highly boosted in the laboratory frame. In our scenario (Table [1]), the lightest 
sbottom has a mass of m^^ = 400 GeV and its boost distribution peaks at ~ 0.95. 

The asymmetries which are based on epsilon products S, see Eqs. ([7]) and (jSj), 
are independent of the sbottom boost, since they are by construction Lorentz in- 
variant. In contrast, the triple product asymmetries are not boost invariant. In 
Figure [6] (left), we show their boost dependence for the three possible triple prod- 
uct combinations T = (^^1^2), (^^1^2), and {Mi). The asymmetries are calculated 
for pp — )■ bib\ production at the LHC, with the subsequent decays 61 — )■ txi and 
Xi — )■ ^i^i. In the sbottom rest frame, = 0, they coincide with the correspond- 
ing epsilon product asymmetries. Note that the size of the asymmetries strongly 
depends on the choice of momenta, and largest values are obtained for the optimal 
triple product T = (^^1^2), as given in Eq. ([9]). 
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Figure 6: Triple product asymmetries A(tiii2), A{biii2), and A(tb£i), see Eq. (|8]) 



for sbottom decay, {bi — )■ txi , followed by Xi 



and t —J- bh'ei2, see Figure [T]), 



as a function of the sbottom boost (left). Asymmetry A{tiii2) in the sbottom 
rest frame and the laboratory frame (at the LHC with ^/s = 14 TeV) as a function 
of (right). The SUSY parameters are given in Table [H 



3.1 Triple product asymmetries in the laboratory frame 

The size of the triple product asymmetry in the laboratory (lab) frame of the LHC, 
is obtained by folding the boost dependent asymmetry A{Pi) with the normalized 
sbottom boost distribution |16] . 

= - APl) df3i, (16) 

of the production cross section a = a{pp — )■ bmb^^). The folded asymmetry A}^"" is 
almost reduced by half, compared to A{[5i = 0). This can be seen in Figure [6] (right), 
where we show both asymmetries as a function of the CP phase (pA^- Thus the 
sbottom boost effectively reduces the asymmetries, but does not change their shape 
with respect to the phase dependence. 

In the following, we will quantify the expected measurability of the the folded 
triple product asymmetry A^^^{tiii2), Eq. (|T6|) . in the lab frame. Since a measure- 
ment of the asymmetries also depends on the production cross section, the sbottom 
boost distribution, and the size of the corresponding branching ratios for sbottoms 
and charginos, we will present the upper bounds on the expected luminosity. 

As defined in Appendix [G], £ is the minimal required luminosity to observe a signal 
at 95% CL above statistical fluctuations. In Eq. (fT7|) . the combined cross section of 
sbottom production and decay is denoted by a. 



11 



In 



nil 










X. 1 






1 


\io\ ; 




1 






■ . . . ■ 


^^^^^ 

. 







600 



500 

> 
O 



400 



300 



mt,<m~x±+mt 



300 400 500 
[GeV] 



500 



400h 
£300 



200 



njl n 3;r/2 2n 
(f>A„ 



look 




m^,<m^f+mt 



^1 



10 



100 



1000 



100 200 300 400 500 
[GeV] 




600 



3.0 

2.5 

^2.0 
> 

^1.5 

.a 

- 1.0 
0.5 
0.0 




10 15 
tan/3 



20 



Figure 7: Contour lines of the minimal required luminosity C, Eq. f[T7j) . in units 
of fb"^ to observe the triple product asymmetry A^^^(tiii2), Eq- f|T6l) . at la above 
fluctuations at the LHC, with y/s^ = 14 TeV, for sbottom production and decay, 

bi txi , followed by Xi ~^ t ^ b Ui £2 ■ The SUSY parameters are given 

in Table [H 



In Figure [71 we show contour lines of the minimal required luminosity C in 
different planes of the MSSM parameter space for the benchmark scenario as given 
in Table [T] The typical sizes of the sbottom and chargino branching ratios BR(6i — >■ 
txf), and BR(x^ — )• ef z/g), are of the order of 15% to 35%. As can be seen, a large 
part of the parameter space can be probed with £ = 10 fb~^. 

Note however that the presented values must be regarded as absolute lower 
bounds on the minimal required luminosities, only. The definition of C, as in 
Eq. f lTT]) . is purely based on the theoretical signal rate and its asymmetry, since 
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detector efficiency effects and contributions from CP-even backgrounds are not in- 
cluded. However, the minimal luminosity requirement can be used to exclude those 
MSSM parameter regions which cannot be probed at the LHC. Clearly, in order 
to give realistic values of the statistical significances and required luminosities, a 
detailed experimental study is necessary, which is however beyond the scope of the 
present work. In particular, the feasibility of the event reconstruction has to be dis- 
cussed. First ideas to attempt the reconstruction of a stop decay chain have recently 
been reported [17] . 

4 Summary and conclusions 

We have analyzed CP observables in the two-body decays of a light sbottom 

h^t + xi- (18) 

The CP-sensitive parts appear only in the top chargino spin-spin correlations, which 
can be probed by the subsequent decays 

X7 + T>h ^ Xi + i^e, ^ = e, /i, 

t ^ b + W, W -^Ui + i2- (19) 

Due to angular momentum conservation, the decay distributions of the final state 
momenta are correlated to each other. Asymmetries of triple products of three 
spatial momenta, as well as epsilon products of four space-time momenta are ideal 
tools to probe the CP-sensitive spin-spin correlations. The CP observables we have 
proposed are sensitive to the CP phases of the trilinear coupling parameter Ab, and 
the higgsino mass parameter n, which might be present in the sbottom and chargino 
sector of the MSSM. 

We have analyzed the asymmetries and event rates in a general MSSM frame- 
work. For maximal sbottom mixing, the asymmetries reach up to 40%. The class 
of asymmetries which are based on triple products are not Lorentz invariant, and 
thus are reduced by a factor of about three when evaluated in the laboratory frame 
at the LHC. Luminosities of at least 10 fb~^ are required to observe a CP signal of 
la above statistical fluctuations at the LHC. 

Clearly, the measurability of the asymmetries and thus of the CP phases can 
only be addressed properly in a detailed experimental analysis, which should take 
into account background processes, detector simulations and event reconstruction 
efficiencies. We want to stress the need for such a thorough analysis, to explore the 
potential of the LHC to probe SUSY CP violation. 

Interrelations between the T-violating electric dipole moments (EDMs) and the 
possible SUSY CP phases underline the need to determine CP observables outside 
the low energy sector, in particular by measurements at the LHC. Since the proposed 
CP asymmetries at colliders depend in a different way on the SUSY parameters 
than the EDMs, they would be the ideal tools for such independent measurements. 
Depending on the experimental results, the EDM bounds could be either verified, 
or the CP-violating sectors of the underlying SUSY models have to be modified. 



13 



Acknowledgments 

We thank F. von der Pahlen for very helpful comments and discussions. This work 
was supported by MICINN project FPA. 2006-05294. The authors thank the Aspen 
Center for Physics for hospitality. 



Appendix 

A Sbottom mixing 

The masses and couplings of the sbottoms follow from their mass matrix [50] 

4/ = -{bh h)\ " I I h (A.i; 



e^'^^mftlA^I m 



R 



)l 









with 



mf^ = M? + (^-^ + ^ sin^ ejj m\ cos(2/3) + m^, (A.2) 

= M\- \n\ sin^ cos(2/3) + m^, (A.3) 

with the soft SUSY-breaking parameters Mg, M^,, the ratio tan /3 = f2/i'i of the 
vacuum expectation values of the two neutral Higgs fields, the weak mixing angle 
the mass of the Z boson, and the mass of the bottom quark. The CP 
phase of the sbottom sector is 

<\>l = arg[As], (A.4) 

A^ = Ab-fi*tan(3, (A.5) 

with the complex trilinear scalar coupling parameter Ab, and the higgsino mass 
parameter Note that for \Ab\ ^ |//| tan/3 we have 0^ ~ (pA^- The sbottom mass 
eigenstates 

are given by the diagonalization matrix 



Q^<Pl cos 01 sin 6^^ 



n'=\ , (A.7) 



sin 9^ e ^'^^ cos 6^ 



with the sbottom mixing angle 9\ 



b 



-mb\A 



2 2 

-I "^i ~ ''^t 

cos% = ^= sin% = ^= ''^ (A.8 



m 
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The mass eigenvalues are 



1 

rrir = - 

bi,2 2 



K\+^6'«)T JK^-"^i^)' + 4m2|A^|2 . (A.9) 



B Chargino mixing 

The complex mass matrix of the charginos is [T] 

V mwv 2cosP fi J 

with the SU{2) gaugino mass parameter M2, and the mass of the W boson my/- 
We optain the chargino masses and the couplings by diagonalizing the chargino 
matrix pQ, 

f/*M^±0 = diag(m^±,m^±), (B.ll) 
with the two independent, unitary diagonalization matrices U and V . 

C Lagrangians and complex couplings 

The interaction Lagrangian for sbottom decay 6^ — ^ ^Xi ^ Kn ~^ '^xt ^ given by [HUI 

-^tl^^ = oHL Pr + Pl) Xt K. + h.c, (C.12) 

with Pl,r = (1 =F 75)72, and the weak coupling constant g = e/ sin 6'^, e > 0. The 
couplings are 

^mi = -"^ml Uil + YbTl^ml Ui2, k^mi = ^ml ^j25 (C.13) 

with the sbottom diagonalization matrix TZ'', Eq. flA.Tp . the chargino diagonalization 
matrices U, V, Eq. (IB. lip , and the Yukawa couplings 

V2mvi/smp v^mvi/cosp 



The interaction Lagrangians for chargino decay Xi ~^ t^vf^ or — > £f z/^, are [T] 



= -gU^.^tPLi^n- 9V:dfPLii>* + h.c., i = e,fi. (C.15) 



D Kinematics and phase space 

For the sbottom decay — we choose a coordinate frame in the labora- 
tory (lab) frame such that the momentum of the sbottom points in the ^-direction 

= (E^, 0, 0, Ip^I), (D.16) 
p'^ = {Et,\pt\smet,0,\pt\coset). (D.17) 
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The decay angle 6t = <{pi, Pt) of the top quark is constrained by sin Of^"^ = \pj^\/\pi\ 
for Ip^I > |p|^| = (m?, m^, m^±)/2mt, with the triangle function A(a, 6, c) = + 
IP' + — 2{ab + ac + be). In this case there are two solutions |39|l58] 



(m? + mf — m^±)|p^| cos6'( ± Ej^./X{m?m'f^ m?±) — 4|p^|2 sin^ 9t 

' ~ 2|p^|2sin2 et + 2ml ' 

(D.18) 

For Ip^I < |p|^| the angle 6t is unbounded, and only the physical solution |p^| is left. 

The momenta of the subsequent decays of the chargino xf ~^ Eq- Q, 

and those of the top quark t — )■ bW, W — )■ i2J^e, Eq. (j4]), can be parametrized by 

p'^ = Ei,{l, sin 9b cos (f>b, sin Ob sirKpb, cos 9b), (D.19) 

= £'£^(1, sin 6*1 cos 01, sin ^1 sin 01, cos 6^1 ), (D.20) 

= £'£2(1, sin cos 02, sin 6^2 sin 02, cos 6^2), (D.21) 

with the energies 



2 2 



-^^1 ~ oTc^ ) \ a — V' ~ 'nTu' \ — \ a — V' (D.22) 

2{E^± - \p^±\ cos9d,) 2{Et - \pt\ cos9dJ 

2 

2{Ew-\Pw\cos9d,)' ^ ■ ' 

The decay angles, 9d^ = <(p^±,p^i), 9d^ = <(pt,Pfe), and 6^02 = <{pw,Pe2), are 

cos6'£,i = p^± ■ P£i, cos9d, =PfPb, cos6'b2 = piy ■ p^2) (D.24) 

with the unit momentum vectors p = p/|p|, and p^± = p^ — Pt, Pw = Pt — Pb- 

The Lorentz invariant phase-space element for the squark decay chain, see Eqs. ([2])- 
(jl]), can be decomposed into two-body phase-space elements [39|[58] 

dUps{si]p^^,Pi^^,Pb,Pe„Pe2) = dUps{si;pt,p^±) 

y<ds^± dUps{s^±]pi^,Pi^^) dst dLips{suPb,Pw) dsw dUps{sw]P£2,Pui), (D.25) 

where we have to sum the two solutions |pf | of the top quark momentum, see 
Eq. (jP.lSp . if the decay angle 9t is constrained. The different factors are 

I ip 12 

dUps{si;pt,p±) = — , ^ ^ — ^sin^ici^t, (D.26) 

Stt \Et \pi\ cos9t - El \pt\\ 

1 \Pi P 

(iLips(s ±;p^,,PpJ = — 3- cil^i, (D.27) 

2(27r)2 — m'^ 

dLips{st;pb,Pw) = tzttt-^ 2 ^ — 2~ ^^b, (D.28) 



2(27r)2 ml — m 



w 



\Pi2\ 



dUps{sw]Pi2,Pu,) = TUTTT^ — i— C^^3, (D.29) 
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with Sj = and dVtj = sinOj dOj d(j)j. We use the narrow width approximation 



m)?ds, = (D.30) 
mjL j 



for the propagators 



A(i) = F' ^^-31) 

Sj — m - + «m,i j 



which is justified for Tj/rrij <^ 1, which holds in our case for particle widths 
Tj < C(l GeV), and masses rrij ^ 0(100 GeV). Note, however, that the naive 
(9(r/m)-expectation of the error can easily receive large off-shell corrections of an 
order of magnitude and more, in particular at threshold, or due to interferences with 
other resonant or non-resonant processes. For a recent discussion of these issues, 
see, for example, Ref. 



E Density matrix formalism 

The amplitude squared for the entire sbottom decay chain, Eqs. ([2])-(jl]), has been cal- 
culated in Ref. [50], by using the spin formalism of Kawasaki, Shirafuji and Tsai [60] . 
We calculate the amplitude squared in the spin-density matrix formalism [5^IM]. 
which allows a separation of the amplitude squared into contributions from spin 
correlations, spin-spin correlations and from the unpolarized part. In that way, the 
CP-sensitive parts of the amplitude squared can be easily separated and identified, 
allowing to find the epsilon product, that yields the largest CP asymmetries. 



In the spin-density matrix formalism of Ref. [52], the amplitude squared of the 
sbottom decay chain, Eqs. can be written as 

|Tp = |A(t)r |A(x-)r \AiW)\'x 

E p^^^ PDAxrf''- PD^f,: pdmkk- (e-32) 

The amplitude squared is composed of the propagators A(j), Eq. (lD.3ip . of particle 
j = ^! oi' and the un-normalized spin density matrices poib), Poiixf)^ PD2{t)^ 
and pD-iiW), with the helicity indices Aj, of the chargino, the helicity indices A(, X[ 
of the top quark, and those of the W boson, A^, A'^. 

The density matrices can be expanded in terms of the Pauli matrices 

H^n''''' i^')x.K ^D, (E.33) 
PoAxtf^'' = S''^'' + {aY^'' J:^^, (E.34) 



PD,{Wf^'^ = Dl'-e'/e^''*, (E.36) 
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with an implicit sum over a, 6 = 1, 2, 3. 

The polarization vectors 5^'= of the W boson fulfil the completeness relation 



Pw,^J, Pw,v 
'9^lv + 2 ' 



m 



(E.37) 



w 



with e^'' = 0. Similarly the spin four-vectors s", a = 1, 2, 3, for the top quark t, 



and s ±, b = 1, 2, 3, for the chargino Xi ? ^Iso fulfil completeness relations 



Ptpt 



m 



2 ' 



and they form an orthonormal set 



/ii^ _|_ A; A-i 

'{J ' 2 ' 

m ± 



(E.38) 



0, (E.39) 



with the notation p'^ = p^/m. 

The expansion coefficients of the matrices, Eqs. flE.33p - flE.35p . are 



D 



y a 



■'D 



9' 



b |2 
mi I 



(+)y 



(+) 2 



2 



+g'Re{lUkl,r} [{st ■ p^){s'^^ ■ Pt) - {st ■ s'^){p^^ ■ p,) 



-/Im{/L(4)*}[4, Pu P^f\ 



yb 



9_ 
2 



IT/ |2/ 2 2 \ 

—\yii\ ("^y± -"^vj, 



Y + - iPb ■ Pt)9n <t)J^e^''''^Pt,a Pb, 



v-ia fif 
D2 



/3 ' 



i^)Y"'t{Hsr+PbSt'' - {Pb ■ s^,)9n ^Ifs^'-'^sl^ p,, 

9^ [p'A+PlP'., - {Pi.-Pu,)9''''] 
with the couplings as defined in Appendix and the short hand notation 

[Pl,P2,P3,P4] = £tMuaP Pi P2 Ps pI; £^0123 = 1- 



(E.40) 
(E.41) 
(E.42) 



(E.43) 

(E.44) 
(E.45) 

(E.46) 

(E.47) 
(E.48) 

(E.49) 
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The coefficients for sbottom decay, Eqs. (lE.40p - (lE.43p . are obtained from those of 
stop decay given in Ref. [16], by the replacements of the couphngs a^j — )■ /^j, and 

— !■ /c^j. The signs in parentheses hold for the charge conjugated processes, that 
is b*^ txt in Eqs. flE^ and flK42D . xt ^ ^t^l in Eq. f lK45|l . t hW in 
Eqs. f lK46|l and ^M), and finally W' in Eq. f lK48|) . 

Inserting the density matrices, Eqs. f lE.33P -f )E.36p into Eq. (IE.32|) . we obtain 

\n = A\m? \mf)? \hw)?x 

[D D, + D, S^f + S^^^ + S^^ S^^ S^f ] Dap., (E.50) 

with an implicit sum over a, b. The amplitude squared |Tp is now composed into an 
unpolarized part (first summand), into the spin correlations of the top (second sum- 
mand), those of the chargino (third summand), and into the spin-spin correlations 
of top and chargino (fourth summand), in Eq. (lE.50p . 

With the completeness relation for the W polarization vectors, Eq. (lE.37p . we 
find 



^D2 -^Spo- 



(E.51) 

(E.52) 



and the sign in parenthesis for the charge conjugated decay, i — )■ bW~ , W~ — )■ u^ll^. 
By also using the completeness relations for the top and chargino spin vectors, 
Eq. flE.38p . the products in Eq. (lE.SOp can be written as 



sr^a sr^a pa 7-) 

^D2 ^3 pa 



It 



{Pb-Pui) 



{P^f ■ Pt){Ph ■ Pt) - m^{p^± ■ pe^) 



lb |2 _ 11,6 



\pab spb 



spa pa p) 
^D2 -'^SpCT 



3_ 
2 



l^.il 



ml±ipt -Ph)- {Py± ■ Pt){Ph ■ Py±) 



9^ ( V^miY + \K, 



\V^l\ 



(E.53) 



(E.54) 



{P^f ■Ph){Pxf ■Pi2)^t 



+ {Pt ■ Ph){Pt ■ Pi2)ml± - {pi, ■ vi^)m\m\ 



-{PX± ■ Pt){Pt ■ Pl2){Pxf 



(Pb-Pui) 



2g'Re{lUkl,r}\Va\'m^±mt{p,-p,^) 

{Pxt ■ Pt)iPei ■PI2)- (Pxt ■ P^2){Pt ■ Pil) 

- 2^® Im{ ( /cl; J * } 1 I ^m^± mt {pb ■ Pu^ )[PhPuPir, PI2 ] , 

(E.55) 

with the short hand notation Eq. (lE.49p . There is no sign change in the terms 
Eqs. (1E.53I) - (]E.55I) for the charge conjugated process b*^ txt^ with the subsequent 
decays xt ^i^h and t bW~ , W~ z/^£^. 
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F Sbottom decay widths and asymmetry 

The partial decay width for the sbottom decay bm txf is [26] 

T{K. ^ txf) = ^ ^^^3 ' D, (F.56) 



with the decay function D given in Eqs. (lE.40p . For the decay we have 

h 

whith the approximation = 0. The sbottom-bottom-neutrahno couphngs are [26] 

^mj — '^ml fbj ~^'^m.2^fji ^mj — '^ml^bj ~^ '^m2 f^j^ (F.58) 

/2 

/5 = ^sin^^(tan^^iV;2-iV;i), (F.60) 

hb, = {KjT = -n(cos/3A73 + sin/3iV;j, (F.61) 

with the sbottom diagonahzation matrix 7^^, Eq. (lA.7p . and the diagonahzation 
matrix for the neutrahno matrix in the photino, zino, higgsino basis, see Ref. [16] 
The sbottom decay width for the complete decay chain, Eqs. ([2])-(jl]), is given by 

T(b ^ vtvtxf^^i^2) = j \T\'^dUps{Sf,;p^^,pj^^,p^±,Pb,Pe„Pe^), (F.62) 

with the phase-space element rfLips given in Eq. flD.25p . 

We obtain an explicit expression for the asymmetry, if we insert the amplitude 
squared |T|2, Eq. f E30|) . into Eq. (P, 

, ^ / Sign(g) S^^ ^ZDs,. dLips 

where we have already used the narrow width approximation of the propagators, see 
Eq. (lD.30p . In the numerator, only the spin-spin terms of the amplitude squared 
remain, since only they contain the epsilon product S, see Eq. ([7]). The other 
terms vanish due to the phase-space integration over the sign of the epsilon product, 
Sign(£^). In the denominator, all spin and spin-spin correlation terms vanish, and 
only the spin-independent parts contribute. Inserting now the explicit expressions 
of the terms Eqs. (1E.40P - (1E.55P into the formula for the asymmetry, Eq. (lF.63p . we 
find 

^ ^ ^ /Sign(g)(pb-p^J [pi,Pt,Pi„Pi2] dLips 
(Pv± -Pei)! {Pt ■ Pi2 ) {Pb ■ Put ) c^Lips 
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with the couphng function 



MC(C.)*} 



1 2 — Tn 2 



/ - - \ ±~"''t 

2 ( Kmil^ + l^miP ) 2m +mt ■^^{^mj(^m,j)*} 
\ / X,- 



(F.65) 



G Theoretical statistical significance 

Assuming that the fluctuations of the signal rate are binomially distributed with 
the selection probability p = 1/2(^+1), the signiflcance of an asymmetry is 



S = -J^=V^, (G.66) 

with the integrated LHC luminosity C, and the cross section a of sbottom production 
and decay as deflned below. The statistical significance S is equal to the number of 
standard deviations that the asymmetry can be statistically determined to be non- 
zero. For example a value of 5 = 1 implies a measurement at the 68% confidence 
level. The minimal required luminosity is then 

C = V4i-lV (G.67) 



The cross section for sbottom production and decay is 

cr = Fat X a{pp hmKn) x Br(6,„ -> txj) x Br(t hW) x Bi{W Uec) 

X Br(xr -> e+z>e) x Br(z>e -> x^e), (G.68) 

for m = 1,2, and i = 1,2. For m, z fixed, the combinatorial factor Fn takes into 
account the possible W and chargino xi decays into leptons with different fiavors. 
We assume that the branching ratios do not depend on the flavor, i.e., Br(iy — )■ 
u^e) = Bt{W — )■ i^^/i), and Br(Xj~ — )■ e'Ue) = Br(x~ — ii'v^). The factor is thus 
Fat = 4, if we sum the lepton flavors e, /x. We further assume Br(i>e — )■ Xi'^e) = 1, 
which applies to our SUSY scenarios considered. 

For the calculation of the sbottom decay widths and branching ratios, we use the 
formulas as given in Eqs. flF.56p . f lF.57p . For the calculation of the chargino decay 
widths and branching ratios, we consider the two-body decays [39] 

xt Be, Ti)^, iiue, fi^2i^r, Wxl, £ = e,/i. (G.69) 

We neglect three-body decays, which are suppressed by phase space. 
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